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Abstract 



In a long wavelength regime, the effective properties of particulate composites, including 
^S| \ nanocomposites, may be estimated using one of various homogenization formalisms, such as 

the Bruggeman and Maxwell Garnett formalisms, and the approach of the strong-property- 
f/3 ' fluctuation theory (SPFT). In the conventional implementations of these formalisms, the con- 

O . stituent particles are treated as point-like scattering centres. However, extended formalisms 

have been established — which involve integral formulations — that take account of the spatial 
extent of the constituent particles. In particular, the extended second-order SPFT takes ac- 
count of both the size of the constituent particles and their statistical distributions. We derive 
explicit representations of the extended second-order SPFT appropriate to isotropic chiral and 
uniaxial dielectric homogenized composite mediums. These results may also be employed in 
^L^. extended versions of the Bruggeman and Maxwell Garnett formalisms. 

Q^, keywords: strong-property-fluctuation theory, depolarization dyadic, Bruggeman formalism, 

Maxwell Garnett formalism 
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'sf ; 1 Introduction 

m 
m 



In the present era of exotic composite materials with nanoscale architectures, there is a pressing need 

^^^ ^ for accurate theoretical tools to predict their electromagnetic properties [1]. While this generally 

Q ■ poses a formidable challenge to theoreticians, matters can be simplified considerably provided 

OO , that wavelengths are sufficiently long relative to the length scales of the particles from which the 

composite is assembled — the composite may then be regarded as being effectively homogeneous. 

Various well established formalisms are available to estimate the constitutive parameters of such 

^ . homogenized composite mediums (HCMs) [2]. One of the most sophisticated — and one which 

H \ has also gained prominence lately in studies of HCMs as metamaterials [3] — is based on the 

strong-property-fluctuation theory (SPFT) [4]. 

The origins of the SPFT lie in wave propagation studies pertaining to continuous random 
mediums [5, 6], but it was later adapted to estimate the constitutive parameters of HCMs [4, 7, 8, 
9, 10]. Unlike other more commonly used approaches to homogenization, such as the Bruggeman 
and Maxwell Garnett formalisms [11], the SPFT can accommodate a comprehensive description 
of the distributional statistics of the constituent particles. The SPFT provides an estimate of the 
HCM constitutive parameters via a recursive scheme, based on an apt ambient medium described 
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by the Bruggeman formalism. Usually the recursive scheme is truncated at the second-order level, 
wherein the spatial distribution statistics of the constituent particles are described in terms of 
a two-point covariance function and its associated correlation length. In a recent development, 
the extended SPFT was established which takes into account the spatial extent of the constituent 
particles [12, 13]. Similarly extended versions of the Maxwell Garnett and Bruggeman formalisms 
have also been developed [14, 15, 16, 17, 18], but these apply only to isotropic HCMs (whereas the 
extended SPFT is available for bianisotropic HCMs [13]) and involve simplistic descriptions of the 
constituent particle distributions. 

A drawback with the SPFT approach to homogenization is that its implementation can be 
an involved process, generally requiring numerical methods. Explicit expressions (i.e., ones not 
expressed in terms of integrals) are not available for SPFT estimates of HCM constitutive param- 
eters, with the exception of certain isotropic dielectric homogenization scenarios [4]. We address 
this issue in the following by deriving explicit expressions for the extended SPFT appropriate to 
isotropic chiral and uniaxial dielectric HCMs. The expressions derived for depolarization dyadics 
can also be utilized in extended versions of the Bruggeman formalism and the Maxwell Garnett 
formalism. 

In the following, vector quantities are underlined. Double underlining and normal (bold) face 
signify a 3x3 (6x6) dyadic. The determinant, inverse and transpose of a dyadic M_ are denoted 
by det [M^], M ~ and M , respectively. The 3x3 (6x6) identity dyadic is represented by / (I); 
and the 3x3 (6x6) null dyadic is represented by (0). Angular frequency is denoted by w; the 
permeability of free space is Ho] and i = ^/—T- The homogeneous bianisotropic medium specified 
by the (frequency-domain) Tellegen constitutive relations [19] 
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is compactly characterized by its 6x6 constitutive dyadic 
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which subsumes the four 3 X 3 constitutive dyadics e, .^ , C. and /i . Subscripts are used to identify 
the particular medium that the constitutive dyadics describe. 

2 Analysis 

Many approaches to homogenization, including those of the SPFT and the Bruggeman and Maxwell 
Garnett formalisms, rely on depolarization dyadics to represent the scattering responses of the 
constituent particles. General integral formulations of depolarization dyadics are presented in §2.1; 
and we show how these are incorporated into the SPFT in §2.2. In §2.3 the main results of this 
communication are presented as explicit representations of the extended SPFT for isotropic chiral 
HCMs and uniaxial dielectric HCMs. 

2.1 Depolarization region 

Let us consider a homogeneous spherical particle of radius rj, embedded in a homogeneous ambient 
medium characterized by the 6x6 constitutive dyadic K . Provided that: 

•' •' ^ amb 



(i) current density distributions induced within the particle are uniform throughout its volume, 
and 

(ii) the particle is small relative to electromagnetic wavelengths, 

the particle's scattering response is captured by the depolarization dyadic [20, 21] 

G (r) d^r. 
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Herein, G (r) is the 6x6 dyadic Green function of the ambient medium. While explicit represen- 
tations of G (r) are not generally available for anisotropic and bianisotropic ambient mediums 
[22], its Fourier transform, namely 



is expressible as [23] 



where 
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Exploiting the spectral representation (4), the depolarization dyadic is given by [20, 23] 
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with q^ = q ' q- 

In order to accommodate a depolarization region of nonzero volume, the depolarization dyadic 
may be considered as the sum [12, 13] 
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and the unit vector q = q/q= {sinO cos<j), sin sine/), cosO). The dyadic D represents the de- 
polarization contribution arising from the vanishingly small spherical region in the limit r/ ^ 0, 



whereas the dyadic D'''(t?) provides the depolarization contribution arising from the spherical re- 
gion of nonzero volume. It is widespread practice in homogenization studies to neglect D"'"(r?) and 
simply take D as the depolarization dyadic [24]. However, the importance of the spatial extent of 
depolarization regions has been underlined in studies of isotropic [25, 26, 14, 16, 17], anisotropic 
[12] and bianisotropic [13] HCMs. 

The integral representation of D has been extensively studied [24]. The volume integral in eq. 
(9) simplifies to the 77-independent surface integral [20, 23] 

D° = — / / G'^Aq) sine dOdc^, (13) 
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which is easily evaluated for isotropic dielectric-magnetic and isotropic chiral mediums [19]. Explicit 
evaluations for uniaxial dielectric mediums have been presented in terms of hyperbolic functions 
[20], whereas an elliptic function representation is available for biaxial dielectric mediums [27]. 

The depolarization contribution given by D^(t?) has been computed in various numerical studies 
[12, 13], but hitherto no explicit evaluations of the volume integral in eq. (10) have been reported 
for ambient mediums other than isotropic dielectric ambient mediums [4]. A notable simplification 
arises in the case of Lorentz-reciprocal ambient mediums (i.e., ambient mediums which satisfy 
e = f"^ ) C = ~C ) s-iid /i = /i"^ [28]). Therein D^(ry) has the surface integral 

— amb —amb = amb = amb ^ amb ^ amb 

representation [10] 
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with K± being the q'^ roots of det 
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and the scalar function 
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Evaluations of the surface integral on the right side of eq. (14) are provided in §2.3.1 and §2.3.2 
for isotropic chiral ambient mediums and uniaxial dielectric ambient mediums, respectively. 



2.2 Homogenization 

Depolarization dyadics play a central role in formalisms employed to estimate the constitutive 
parameters of HCMs [4]. Here we focus on the approach provided by the extended SPFT [12, 13]. 
We concentrate on the homogenization of a two-phase composite wherein the two constituent 
phases, labelled as a and b, comprise spherical particles of average radius 77. The constituent phase 
a occupies the region V^ whereas constituent phase b occupies the region V^,, The constituent phases 
are randomly mixed with their distributional statistics being described in terms of moments of the 
characteristic functions 
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The volume fraction of phase ^, namely /^ , is given by the first statistical moment of ^i ; i.e., 
(<l>^(r)) = fi . We have (^afc)) + (*^b(zi)) = 1- For the second statistical moment of $£, the 
physically-motivated covariance form [29] 
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is commonly implemented. Herein, the correlation length L is required to be much smaller than the 
electromagnetic wavelengths, but much larger than the constituent particle radius r]. Parentheti- 
cally, the specific form of the covariance function has only a secondary influence on SPFT estimates 
of HCM constitutive parameters, across a range of physically-plausible covariance functions [30]. 

The nth-order SPFT estimate of the HCM constitutive dyadic, namely K}^' , derives from the 
recursive refinement of an ambient homogeneous medium, characterized by the constitutive dyadic 
K . At lowest order (i.e., zeroth and first order), the SPFT estimate of the HCM constitutive 
dyadic is identical to that of the ambient medium [10]; i.e., 

kM =kW =K . (18) 

Furthermore, K is delivered by solving the nonlinear equations 
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wherein the polar izability density dyadics 
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and K „ (^ = a,b) are the constitutive dyadics of constituent phases a and b. 

For a broad range of bianisotropic HCMs, the SPFT (including the extended SPFT [31]) con- 
verges at the second-order leveP [32] . The second-order SPFT estimate of the HCM constitutive 
dyadic is [10] 
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with the mass operator dyadic term [5] — corresponding to the covariance function (17) — being 



S[2](r?,L) = /a/Jx {v)-X (r?)l •D+(L)- \x (v) - X iv) 
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Thus, within the extended second-order SPFT, the estimate of the HCM constitutive dyadic de- 
pends on two length scales: the constituent particle size rj via the depolarization dyadic and the 
correlation length L via the mass operator. 

2.3 Explicit representations 

We now turn to the evaluation of the surface integral which delivers D"*". This term crops up as 
D'''(t?) in the extended depolarization dyadic representation (8) and as D"'"(L) in the mass operator 
term (22) which yields the second-order SPFT estimate of the HCM constitutive dyadic. A key 



The second-order SPFT approximation is also known as the bilocal approximation. 



step in the evaluation of the integral on the right side of eq. (14) is the approximation of the 
g«w (^i _ ifjq"j term in the integrand by its asymptotic expansion 1 + (rjq)'^ /2 + i{r]q)^ /3, which is 
permissible since ri^fK± <C 1 in the long wavelength regime. The evaluation of D"'"(L) is isomorphic 
to that of D"'"(7?) since here we similarly have L^/k± ^ 1 in the long wavelength regime. In the 
following we consider two types of ambient medium: an isotropic chiral ambient medium in §2.3.1 
and a uniaxial dielectric ambient medium in S2.3.2. 



2.3.1 Isotropic chiral ambient medium 

Suppose that the ambient medium is an isotropic chiral medium. Its constitutive dyadic has the 
form [33] 
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The integral on the right side of eq. (14) may then be evaluated to give 
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after some straightforward manipulations. 
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2.3.2 Uniaxial dielectric ambient medium 

Suppose that the ambient medium is a uniaxial dielectric medium. Without loss of generality, let 
us take the distinguished axis of this uniaxial medium to be aligned with the Cartesian z axis. 
Accordingly, its constitutive dyadic takes the form [19] 
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Within the present context of the SPFT, such an ambient medium arises in homogenization scenar- 
ios in which one of the constituent phases is a uniaxial dielectric medium and the other is either an 
isotropic dielectric medium or a uniaxial dielectric medium. Hence, we observe from eqs. (20) and 
(22) that only the upper three diagonal entries of D^, namely [P"*"] (n = 1, 2, 3), contribute to 
the estimate of the HCM constitutive dyadic. After some manipulations, these entries are evaluated 



as 
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with the dimensionless scalar 



(26) 



(27) 



(28) 



The representations (26) and (27) apply when 'y^„,\, is complex-valued (with nonzero imaginary 
part). If 7amb is real-valued then the representations (26) and (27) apply when < 7amb < 1; 
for 7amb < 0, the ^/y^ tanh" (■y/7amb) term in eqs. (26) and (27) should be replaced by 
— ■^— 7amb tan^-*^ (■^— 7amb)- In the scenario 7amb > 1 — which corresponds to nondissipative uniax- 
ial dielectric ambient mediums with indefinite permittivity dyadics [34, 35, 36] — the components 
of D~''(r/) are undefined, as are the corresponding components of D [20]. 



3 Discussion 

The main results of this communication are the derivations of the eqs. (24), (26) and (27) which, 
when substituted into eq. (21), yield explicit formulations of the extended second-order SPFT for 
isotropic chiral HCMs and uniaxial dielectric HCMs (when supplemented with the corresponding 
expressions for D which are available elsewhere [19]). For more complex HCMs, numerical methods 
are needed to evaluate the depolarization dyadic D(ry) and mass operator S[^l(r/, L). 

It is helpful to consider the isotropic dielectric specialization of the extended second-order 
SPFT result (21). In this case, the constituent phases are both isotropic dielectric mediums with 
permittivities e^ and e^^; in consonance, the ambient medium is also an isotropic dielectric medium 
with permittivity e^^y^- The corresponding second-order SPFT estimate of the HCM permittivity 
is 

J2] _^ 1 ( Sl^](r?,L) ^ 
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wherein the depolarization scalar may be expressed as the sum 
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The contribution to d{r]) associated with a vanishingly small depolarization region is provided by 
the well known result [19] 

d' = ^:r^—, (31) 



whereas the contribution associated with a depolarization region of nonzero volume may be ex- 
tracted from eqs. (26) and (27) in the limit e 
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amb 
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This expression for df^ij]) also follows from eq. (24) when C,^^^ = and /iamb = Mo- The scalar mass 
operator term in eq. (29) is provided by the isotropic dielectric specialization of eq. (22) as 

S[21(7?,L) = /J,[Xa(^)-Xb(r/)]'d+(i.), (33) 

with the polarizability density scalars 
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We note that the explicit expression for the extended second-order SPFT estimate of the HCM 
permittivity (29) is consistent with a corresponding result derived for the unextended SPFT, within 
the long wavelength limit represented by [4] 

^''1 -^ --SP](L), (35) 
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where SP](L) = SP](?7,L) evaluated with d^ in lieu of d{r]). 

While our focus has primarily been on the extended second-order SPFT, we bear in mind that 
the expressions (24), (26) and (27) can also be deployed in extended versions of other homogeniza- 
tion formalisms, such as the frequently used Bruggeman and Maxwell Garnett formalisms. In fact, 
the Bruggeman homogenization formalism is equivalent to the zeroth order SPFT [9]. Hence the 
extended Bruggeman estimate of the HCM constitutive dyadic is simply the dyadic K which 
may be extracted from eq. (19). Like the Bruggeman formalism, the Maxwell Garnett formalism 
is based on depolarization dyadics [24], but with one of the constituent phases playing the role of 
the ambient medium. Therefore, the Maxwell Garnett formalism — including its incremental [37] 
and differential [38] variants — may be extended by implementing the appropriate ?7-dependent 
depolarization dyadic Dfr?). 

We close by considering the question: For what range of constituent particle sizes are the 
extended homogenization formalisms discussed herein appropriate? An upper bound is straight- 
forwardly established by the requirement that the particles must be small relative to the electro- 
magnetic wavelength(s), in order to be consistent with the notion of homogenization. For optical 
wavelengths, the linear dimensions of the constituent particles must therefore be at most 38-78 
nm. A lower bound comes into effect because, at sufficiently small length scales, quantum pro- 
cesses cannot be neglected in the description of the constituent particles and their interactions. 
Accordingly, the lower bound is material-dependent. We note that for very small constituent par- 
ticles, their constitutive parameters may differ significantly from those of the corresponding bulk 
materials and depend upon the shape and size of the particles [39]. This is particularly the case 
for metallic particles smaller than the mean free path of conduction electrons in the bulk metal, 
wherein the mean free path may be dominated by collisions at the particle boundary [40]. Within 
this particle-size regime, a recent study using spectroscopic ellipsometry demonstrated that an ex- 
tended Maxwell Garnett homogenization formalism adequately characterizes a HCM based on silver 
nanoparticles as small as 2.3 nm [41]. This study also highlights the prospects of implementing 
extended homogenization formalisms, such as those described herein, in particle sizing applications 
for nano composites. 
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